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Abstract 



We examine angular distributions of leptons arising from the decay 
of J/ip in inclusive hadroproduction. Taking into account feed-down 
contributions from xi,2, a flat distribution emerges without fine tun- 
ing parameters. Changes in the ratio of direct to total J/ijj cross 
. sections would change this distribution significantly. Such angular 

distributions are, therefore, confirmatory tests of the J/ip suppression 
signal for the production of a quark-gluon plasma. Related effects are 
predicted for the T. 
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1 Introduction 



In the hadroproduction of J/i/j, the angular distribution of the dileptons 
coming from the decay of the quarkonium is of the form 

l + acos 2 6 (1.1) 



dcosO 



where 9 is the angle (measured in the rest frame of the quarkonium) between 
the direction of motion of the positively charged lepton and the (lab-frame) 
direction of the quarkonium momentum. As long as the initial particles 
are unpolarized, there can be no difference in the weights for production of 
quarkonia with helicities m and —m. The angular distributions can then be 
summarized by an "alignment parameter" Jl| — 

C = — , and a = r- (1-2) 

<j l + c 

Here <tl is the cross section for production of quarkonia with longitudinal 
polarization (J z = 0), and a is the total cross section. For hadroproduction, 
a is close to zero for all beam and target combinations investigated P, §], 
implying that ( is close to 1/3. 

The question of angular distributions and alignment has attracted a lot of 
attention recently. In ||] it was pointed out that heavy-quark spin symmetry 
has strong implications for (. Non-relativistic QCD (NRQCD) ||, which 
incorporates heavy-quark spin symmetry, is a framework for understanding 
production rates of quarkonia. Polarization in hadroproduction has been 
investigated in NRQCD in || [7|, |^]. The dependence of a on kinematics is 
claimed to be a discriminant between different models for hadroproduction of 
quarkonia, and hence a crucial test of NRQCD |J. Polarization of quarkonia 
produced in the decay of Z has also been investigated fllP] , and may be 
tested at LEP. 

The phenomenological application of NRQCD to inclusive hadroproduc- 
tion of J /ip has been quite successful. Absolute cross sections are fully com- 
patible with the scaling laws of NRQCD, and the energy dependence of total 



cross sections g, [LI]], as well as longitudinal momentum distributions |L2 
are reproduced rather well. The ratio of Xi an d X2 cross sections can also be 
understood by going to higher than leading orders in the NRQCD expansion 
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parameter, v [13|]. In view of this, the alignment of quarkonia produced in in- 
clusive hadron-hadron collisions is investigated. Although the computations 
are performed to leading order in v 2 for each channel, this already leads to 
an expansion to order v 9 in some cases. 

Our results are summarized here. The observation of almost vanishing a 
in inclusive J / ip production arises from the mixture of direct J / ip and feed- 
down from x states. We compute alignments in NRQCD at sufficiently high 
order, thus supplementing the results of [0, § by formulae for ( x . The value of 
a depends only on the ratios of certain non-perturbative matrix elements, and 
is well constrained by the scaling laws of NRQCD. A flat angular distribution 
emerges because of the feed-down from x states. 

This has implications for the J/ip suppression signal of the quark-gluon 
plasma [14]]. A plasma melts not only the J/if), the more massive quarko- 



nium states are disrupted even more easily — indeed it has been speculated 
that some of them dissolve even at temperatures lower than the QCD phase 
transition Hl5| . If this is so, then a suppression of the J/ip cross section due 



to a plasma must be coupled with a large change in a. This provides a confir- 
matory test for the J / ip suppression signal of the formation of a quark-gluon 
plasma in relativistic heavy-ion collisions. 

In section || the computation of the angular distribution in NRQCD is 
outlined in brief. Details of the computation are given in appendix |A], where 
a helicity technique used in conjunction with previously known methods of 
computation is documented. In section || the results for ( are summarized. 
In the final section [| these results are applied to inclusive hadroproduction 
and the signal of quark-gluon plasma formation. 



2 Non-relativistic QCD 

NRQCD is a low-energy effective theory for quarkonia [Q. The action is 
written in terms of all possible operators consistent with the symmetries of 
QCD. All momenta in NRQCD are cut off by some scale A, taken to be of 
the order of the heavy quark mass, m. The coupling associated with each 
term is identified through a perturbative matching procedure. In NRQCD 
the perturbative short-distance (< 1/A) physics of the production of a heavy 
quark pair (QQ) is factored from the non-perturbative long-distance physics 
of its hadronisation to a quarkonium. A proof of factorisation has been given 
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for production at large transverse momenta and successful phenomenology 



has been done [16|. An explicit test of factorisation for inclusive cross sections 



is provided by a recent next-to-leading order, a§, computation JT7| . 

The NRQCD factorisation formula for the inclusive production of heavy 
quarkonium resonances H with 4-momentum P can be written as 

d " = 5(S^i:c«(K,n(/f)K3), (2.i) 

where $ is a flux factor. The coefficient functions are computable in 
perturbative QCD and hence have an expansion in the strong coupling a s /7r 
(evaluated at A). Although each matrix element in the sum above is non- 
perturbative, it has a fixed scaling dimension in the quark velocity v. Then 
the NRQCD cross section is a double series in a s /ir and v 2 . For charmonium 
a numerical coincidence, a s /ix < v 2 , makes it necessary to consider higher 
orders in v 2 before going to higher orders in a s . 

The fermion bilinear operators /Q are built out of heavy quark fields 
sandwiching colour and spin matrices and the covariant derivative D. The 
composite labels % and j include the colour index a, the spin quantum number 
S, the number of derivative operators N, the orbital angular momentum L, 
the total angular momentum J, the helicity J z and the parity. The hadron 
projection operator 

11(F) = "£\H,s){H,s\, (2.2) 

s 

(where s denotes hadron states with energy less than the NRQCD cutoff), 
is diagonal in appropriate bases j|, ||, 0. When the final state helicities 
are summed it is diagonal in J, J z and parity. When the hadron helicity is 
observed, heavy quark spin symmetry makes it diagonal (up to corrections 
of order v 2 ) in the basis {S, L, J z }. 

The J z -dependence of these matrix elements can be factored out using 
the Wigner-Eckart theorem — 

</Qn(#)/cl> = i- 2 o»( 2S ^L»), 

J 1 (2-3) 
</CTI(tf)/Cj> + h.c. = jjV^^L^^Lf,). 



The factors of J = \/2J + 1 come from Clebsch-Gordan coefficients and are 
conventionally included in the coefficient function. NRQCD supports a power 
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counting rule for each matrix element — 

{K$l{H)K)) = R H A D v d , (2.4) 

where D is the mass dimension of the operator, and the NRQCD order, d, is 
given by the rule 

d = 3 + N + N' + 2(E d + 2M d ). (2.5) 

Ed and Md are the number of colour electric and magnetic transitions re- 
quired to connect the hadronic state to the state /Q|0). At tree level the 
dimensionless number Rh in eq. ( |2.4[ ) must depend only on the hadron H 
under consideration. At higher orders in a s , Rh may be corrected by loga- 
rithms of A. 



Matrix Elements 


Values 


Rh 


Of 4 ' ( 3 5 x ) 
O 3 ^ (%) + ^0 J f ( 3 P ) 


(1.12 ±0.14) x 1(T 2 GeV 3 
(3.90 ± 1.14) x 10~ 2 GeV 3 


0.22 ±0.03 
0.17±0.05 


Of ( 3 S X ) 
Of (%) + ^Of ( 3 P ) 


(0.46 ± 0.08) x 10~ 2 GeV 3 
(1.60 ±0.51) x 10~ 2 GeV 3 


0.09 ±0.02 
0.07 ±0.02 



Table V. The values of various non-perturbative matrix elements in the non- 
relativistic normalisation fllS] . These can be converted to the relativistic 
normalisation by multiplying by 4m. The dimensionless number Rh (eq. 
has been extracted using A = m = 1.5 GeV and v 2 = 0.3. 



In principle, we would like to do phenomenology with experimental mea- 
surements of the matrix elements needed. However, the large number of 
different matrix elements required for fixed target phenomena makes this 
impossible at present. Consequently, we are forced to use the scaling laws 
of eq. ( |2.4j ). Any number obtained in this way can at best be indicative, 
and has to be justified by more detailed numerical work when better data 
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becomes available. The first step is, of course, to test eq. ( |2.4| ). From the 
data summarized in Table [l| it seems that tree-level NRQCD scaling can be 
accepted working hypothesis. 

The coefficient functions, Cy, are computed using a non-relativistic de- 
composition of heavy quark spinors and a Taylor expansion of the matrix 
elements, Ai, in the relative momentum, q, of the QQ pair ||. Spherical 
tensor techniques are then used to recouple the 2-component spinor bilinear 
operators fll3| . The projection to specific J z components can be performed 



at the matrix element level by appropriate choice of gauge. The process of 
matching the perturbatively computed \A4\ 2 to the NRQCD formula of eq. 
( |2.1|) is then simple. 



3 Quarkonium alignment 

We need to compute the alignment parameter only for J/ip, Xi an d X2- 
For completeness, we also list it for the yet unestablished 1 P\ state h c . A 
simple extension of the results in || is sufficient to show that = for all 
quarkonia H to all orders in v. The total cross section has been listed in JTj| 
to order a 2 s v 9 . In this section only the longitudinal parts of the gg — > QQ 
cross sections are listed. Some details are given in Appendix [A]. 

3.1 Direct J/i/j alignment 

The direct J j ip subprocess longitudinal cross section is 



ll4^ (7) + {l^^ (9) + ^ (9) } 



(3.1) 



where 



ir 3 a 2 



lf= 4^ (J_4m) ' (3 ' 2) 
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and 9^^{n) denotes combinations of non-perturbative matrix elements from 
the colour amplitude a (= S, D or F) at order v n . These can be written as 



5m 6 V 3 

n 

2m 6 

J/i> 



(3-3) 



The 9jj f (7) term agrees with that given in ||. The rest of the terms did 
not appear there because the Taylor expansion was truncated at lower order. 
Terms involving l P\ operators, given in M, contribute at order v 11 and hence 
are not included in our results. 

The leading, order v 7 , terms have been considered before || [/J. Under 
the assumption of tree-level NRQCD scaling this gives 

CO* = l + 0(v 2 ). (3.4) 

Order v 2 corrections come from two sources — the terms in the production 
cross section given above, as well as from corrections to the NRQCD action 
through terms that violate heavy-quark spin symmetry. 

3.2 xi alignment 

The longitudinal cross section for xi is 



99 ' 

with the combinations 



L ^(9) + ^(9) + ^(9) 



(3.5) 



^(9) = ^4<Wo) 

02(9) = ^f^SD + ^On^) ( 3 - 6 ) 

^(9) = ^o?CPl) + ^o?( 3 sf) + -^or( 3 D!). 

bm* 18m° 18m° 
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Although Xi production begins at order v 9 , the large number of terms in- 
volved makes its cross section comparable to that of Xo,2- Alignment of Xi 
has not been considered before in the literature. Under the assumption of 
tree-level NRQCD scaling we have 



3.3 X2 alignment 



The leading term in the total cross section for X2 production is given by a 
J z = ±2 operator scaling as v 5 . The most significant J z = operator scales 



as v 7 - 



1 . 5 



„,„,,, ^0f(7), and 9f(7) = -^p?(*Pl*I#) (3.8) 

Since the total cross section starts at order v 5 , the alignment is 

Qg = + O(v 2 ). (3.9) 

3.4 h c alignment 

The production cross section for the 1 P\ charmonium state is— 

5 



144 



0&(5)+6fc(7) 



(3.10) 



where the combinations of non-perturbative matrix elements can be written 

as 

1 23 

^ (5) = 2^°^^ and ^ (7) = "^3W^ C(l5 °°' l5 ° 2) - 

(3.11) 

In this case, 

Cgg = \ + 0(V% (3.12) 
implying that angular distributions are trivial. 
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4 Phenomenology 



In computing the effective alignment of J/ip, the feed-down from radiative 
decays xj iJ/ty must be taken into account. The total spin-projected 
cross section for J/ip can be written as 



(4.1) 



M 



where Bj denotes the branching ratio for the decay of xj to J/ip, and Bjpj^ M 
is the branching fraction for a xj with spin projection M to give a J/ip with 
spin projection J z . We use the notation 



fj = 



B J a x., 



&j/4> + Bi<J Xl + B 2 <7 X2 ' 
to write the effective alignment parameter as 



(4.2) 



c = 0^(1-/1-/2)+/! 



PooCxi+^oia-CxJ 



+ /2-P0V 



(4.3) 



±2 



Po,j+Po,~j 



In writing the equation above, the fact that X2 is produced only with J ; 
at leading order in v is taken into account, and the notation Pq j 
is introduced. Observe that the first term on the right has the form of a 
dilution factor over (j/^, driving the effective ( away from the required value 
of 1/3. This has to be compensated by effects of the other two terms. 

The cross section for the production of any charmonium state is ob- 
tained by convoluting the parton level cross sections with parton luminosity 
factors — 



<J H {J z ;^fS) = C, q (^)a»(J z )+£ gg (VS)a»(J z ), 



(4.4) 



where v^S is the centre of mass energy at which the cross section is measured. 
Recalling that the alignment from the qq — > QQ channels is identically zero, 
it is clear that 



~ Cgg 



1 _|_ ^ m a<iq 

t-gg u gg 



(4.5) 



The second factor is a dilution. The ratio of parton level cross sections on the 
right is determined by the NRQCD scaling laws. As a result, this dilution 
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Figure 1: Measurements of B 2 fi/Bif 2 (collected in [[J) in up (filled circles) 
and pp collisions (open circles) are compared with the predictions of NRQCD. 
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factor is completely determined once the parton densities are specified. In 
the rest of this paper the GRV LO parton densities are used for both proton 
and pion [T9| . 



The quantities fi and /2 have been investigated in several experiments. 
Pion beams have been used at many energies, and there is weak evidence 
for a non-trivial y/S dependence. The number of measurements with proton 
beams is more limited. Some of the older experiments, at low y/~S, ]20 had 



reported /1//2 somewhat smaller than those observed with pion beams. A 
new measurement at higher y/S pTI seems to give a value consistent with 



the pion data at a similar energy. The data has been collected in ||. Using 
the cross sections for \ production computed to order a s v 9 |L3|], and tree 



level NRQCD scaling (eq. |2.4p with R Xl = R X2} reasonable agreement with 
data is obtained^. As shown in Figure [I], NRQCD scaling implies that the 
/1//2 ratio should show significant beam dependence at small yfS but should 
rapidly converge to a smaller common value with increasing yfS. This is due 
to the relatively large influence of Cq q at small y/S, specially for 7r-A collisions. 
Since the main uncertainty is due to a s corrections to eq. (|2.4|), we estimate 
a theoretical uncertainty of about 25% for predictions of fi/ }2- 

The next step is to constrain the three parameters pg , P \ and P 2 2 . In the 
decay x.j ~^ jJ/^i the photon energy is of the order of 0.5 GeV, and much 
less than the NRQCD cutoff. As a result, the decay matrix elements, and 
hence p$ M , cannot be computed in NRQCD through the usual factorisation 
approach. One way out is to obtain them by direct measurement. Each of the 
LEP experiments should have a few hundred identified x' s m the hadronic 
decay mode of the Z . The parameters p$ M can be estimated from angular 
correlations between the decay products in the chain \ —>■ jJ/ip — * 7e + e~. 

An alternative technique has been used in the literature [p2]j . Assuming 
that the decay is dominated by an electric dipole transition, using heavy- 
quark spin symmetry, and integrating over the direction of the photon mo- 
mentum, the parameters can be shown to be just squares of certain Clebsch- 
Gordan coefficients — 

pLm = \C(1J1;M-J Z ,M)\ 2 . (4.6) 



2 For some colour singlet matrix elements, Rh is related to the radial wavefunction 
of the hadron H . Although this intuitive picture is violated in NRQCD, in particular 
by logarithms of A, it provides the motivation for the choice R Xl = R X2 . It would be 
preferable to use data, when it becomes available, to hx these two parameters. 
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Contributions of higher multipoles are expected to be subdominant, since 
they go as powers of E 7 /Mj/^ (E 1 is the energy of the decay photon, and 
Mj/^ is the mass of the J/ijj). With this model, pl = 0, P \ = 1, Pq 2 = 0. 
In it is reported that in 7r-Be collisions at y/S = 31.8 GeV, 

a = -0.01 ±0.12, (4.7) 

implying £ = 0.34 ± 0.05. The same experimental collaboration also reports 

B f 

f l + f 2 = 0.443 ± 0.041 ± 0.035 and = 0.57 ±0.18 ±0.06. (4.8) 

Extracting the values of fi and f'2 from these measurements, we find that eq. 
(|4.3| ) and the assumption of dipole dominated decay give 

C = 0.30 ±0.04 and a = 0.08 ± 0.08. (4.9) 
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The errors in the prediction are obtained solely by propagation from errors 
in fi and f 2 (statistical and systematic errors are added in quadrature). We 
have assumed that the errors in fx and f 2 are uncorrelated. The theoretical 
uncertainty from neglecting order v 2 corrections to ( gg and order a s correc- 
tions to tree level NRQCD scaling is not displayed in eq. ( |4.9|) . 

We can also investigate the allowed departure from the dipole model by 
taking f\, f 2 and a from the above experiment, and finding the contours of 
allowed pg , P \ and Pg 2 using eq. ( |4.3| ). The 1-a allowed contours are shown 
in Figure || It is clear that the dipole model, as well as large deviations from 
it, are allowed. 

Similar results follow if we combine measurements of a and f± j2 from 
different experiments with the same beam and roughly similar energy. We 
can also make these estimates by combining the NRQCD scaling predictions 
for /1//2 with experimental measurements of fi + f 2 - The same features 
emerge — the dipole dominance model of the radiative decay of x used in 
conjuction with eq. (|4.3| ) works well, and departures from it are allowed. In 
other words, strong tuning of parameters is not required to reproduce the 
angular distribution of leptons from the decay of hadroproduced J ftp. 

An application to bb systems may also be considered. In this case it is 
estimated that v 2 ~ 0.1. As a result, hadroproduction of T should have 
roughly equal contributions from the order v 5 production of \2 (with B 2 = 
0.2) and the order v 7 direct production of T. The other \b states should be 
much less important — the xo contribution being suppressed by the smaller 
branching ratio, B < 0.06, and the Xi by the relative order v 4 in the matrix 
elements. Assuming therefore, f 2 = 0.5 and f\ = 0, eq. (|4.3| ) predicts «x ~ 
0.25-0.70. The larger value is preferred if Xb decays occur through electric 
dipole transitions. 

Consider now nucleus-nucleus collisions where the central rapidity region 
reaches a temperature sufficient to suppress x but not J/ip Jl5|, i. e., fi(T) « 
f 2 (T) « 0. The total J/ip cross section should decrease by about 50%. 
Simultaneously the angular distribution should change, with a new ( = 3/7. 
As shown in Figure the change should be easily visible. 

Thermal effects on T have not been considered earlier, since this resonance 
is not expected to dissolve at the temperatures which the LHC may reach 



W\. However, the Xb states are expected to disappear at such temperatures 



15(1 - O ur earlier analysis would then lead us to believe that some reduction 



of the T cross section may occur, and the angular distribution may change to 
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Figure 3: The parameter a as a function of the fraction of J/ip produced by 
decay of xi,2- The dark band is the prediction for f = fi = fi- The lighter 
horizontal band is the 1-a range for a seen in T = experiments. 
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give «x = 3/7. Since the successive disappearance of different quarkonium 
states can be used for thermometry of the plasma, we believe it is important 
to look for this effect. 

We present a list of experiments which might be used to constrain and 
test the applicability of NRQCD to inclusive quarkonium hadropro duct ion. 

• More refined measurements of the cross sections for xi an d X2 pro- 
duction with proton and pion beams, at y/S « 40 GeV, should test 
the scaling of the matrix elements and fix the ratio Rj/^/R x . High- 
statistics experiments at y/S < 25 GeV designed to see the difference 
in /1//2 for proton and pion would be welcome. 

• Measurement of angular correlations in the decay \ — > 7e + e~ would 
test the dipole dominance model of this decay, and hence the NRQCD 
explanation of the angular distribution. Hadroproduction experiments 
are not needed for this. In fact e + e~ machines, such as the LEP, would 
be preferable. 

• Present day errors on a are large, and need to be reduced significantly. 
Simultaneous measurements of /1, fi and a by one experiment can test 
NRQCD scaling rather accurately. 

• The bb system remains almost completely unexplored. Due to the small- 
ness of v 2 in this system, one would expect the phenomenology to be 
quite different from the cc system, and scaling arguments to work bet- 
ter. In view of the importance of T and \b for tests of NRQCD scaling, 
and possible applications to the thermometry of a quark-gluon plasma, 
a vigorous experimental effort in this field would pay large dividends. 

The results can be summarized as follows — an NRQCD based argument 
yields aj/^, ~ 0, provided the scaling formula of eq. ( |2.4| ) is used. The effec- 
tive angular distribution (eq. |4.3| ) needs three parameters related to angular 
correlations in the decays xj ~~ > iJ/i*- They need not be fine tuned to repro- 
duce the observed angular distributions, although the assumption that the 
decays proceed through an electric dipole transition is supported. A crude es- 
timate for the bb system gives ay ~ 0.25-0.70. These findings may be applied 
to the J I ip suppression signal for the formation of a quark-gluon plasma — a 
suppression should be coupled with a large change in a, since the feed-down 
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processes disappear rapidly with increasing temperature. Suppression of Xb 
should similarly influence the angular distribution of T. 
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A Computing the coefficient functions 



We take the momenta of the initial particles, p\ and P2, and the momentum 
of the pair, P = pi + P2 {P 2 = M 2 ), to lie along the ^-direction. We choose 
p\ > and take this to be the axis of quantization of angular momenta. The 
4- momenta of Q and Q {p and p respectively, p 2 = p 2 = m 2 ) are written as 

p = ^P + Lrf and p = \p - Lrf . (A.l) 

The space-like vector q is defined in the rest frame of the pair, and boosts 
it to any frame. For gg initial states we choose polarization vectors 

e X (A) = -^=(0,-A,-i,0) and e 2 (A) = -L(0, A, -i, 0). (A.2) 

Since the polarization vectors are orthogonal to both the initial momenta, 
this corresponds to a choice of planar gauge. 

Euclidean 3-tensors are expressed as spherical tensors. As an example, a 
3- vector etj is written as a spherical tensor of rank 1, with components 

a±i = =F-^= (a x ± ia y ) , a = a z . (A. 3) 

In this representation 

e^A)-^ = -6 m , x , (A.4) 
where m is a helicity index. An useful identity is 

ajbj = a Q b - (a +1 b_i + a_i6+i) = -V3[a, 6]°, (A. 5) 

We have introduced the notation [a, b\ J M to denote two vectors a and b coupled 
to total rank J and helicity M. The coefficient of the terms [a, 6]q can be 
obtained from the appropriate Clebsch-Gordan coefficients. 

Since we work to lowest order in the QCD coupling, the perturbative 
projector is 

n(QQ) = \QQ)(QQ\. (A.6) 
Its normalisation is completely fixed by the relativistic normalisation of states 

(Q(P,OQ(P,V)\Q(P',OQ(P',V')} = 4E p E p (2ir) 6 5 3 (p-p')5 3 (p-p'), (A.7) 
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where the spinor normalisations are = r\ '77 = 1. Expanding E p = Ep = 
y/m 2 + q 2 in q 2 allows us to write any spinor bilinear in terms of transition 
operators built out of heavy quark fields. 

The matrix element for the subprocess qq — > QQ has been treated to 
leading order in v in ||. The result, (g q = for all quarkonium states, can 
be trivially generalized to all orders in v. The squared matrix element can 
be written in the form 

\M\ 2 = [k 1 -L i k 2 -L j + k 2 -L i k l -L j -k 1 -k 2 L i -L j \W i \ (A.8) 

where W 1 ^ represents the square of a heavy-quark spinor bilinear. In the 
kinematics appropriate for this problem, the prefactor to W %: > becomes (1 — 
SiaSjs). Transforming to spherical tensor components, this can be seen to 
give a vanishing longitudinal cross section. 

The computations for the gg process can easily become tedious. We 
introduce here a helicity technique based on a decomposition of the matrix 
element into gauge invariant amplitudes — 

M = \g 2 5 ab S + \g 2 d abc D c + l -g 2 f abc F c . (A.9) 

These amplitudes have been enumerated earlier JT3[ . 

The invariant amplitudes are most compactly written in terms of two 
quantities A and S. When the initial gluon helicities are A and A', with the 
gauge choice given earlier, we get 

A = -i\dx >X i and S^V = -V3 [a, b)° 5 A , A , + V2 [a, bf 2X 6 X ^ (A.10) 

In order to identify all terms to order v 9 we need the colour amplitude S to 
order q 5 — 

S = - («?p) A (ety + ^S jm (q^rj) - Az m z n (q m q^ V ) 



1 16 - . >' \ c 

\_ i™ 1 n v M + 2m i m n p\ ^ Q Q <, 17 V) 



M 



Az m z n z p z r (q m q n qPq r ^ V ) 



1 64 
+ M 5 



'-'jm^n^p jJJ _|_ QfYi 3 m n P 



z r z s (q m q n q p q r q s ^a j rj). 

(A.11) 
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The amplitude D differs only through having colour octet matrix elements 
in place of the colour singlet ones shown above. For the colour amplitude F 
we need the expansion 

FC = ~ (^?) (q m ?T c V ) + ^S jm z n (g'VeVT^). (A. 12) 

In all three colour amplitudes, the terms in A are spin singlet and those in 
S are spin triplet. 

The structure is very simple. Conservation of angular momentum shows 
that the J z = terms are obtained when A = A'. With eq. (|A.10|) recoupling 



of all spherical tensors into terms of well-defined L, S and J can now be 
performed at the amplitude level. This also simplifies the computations 



presented in flTj|. The full procedure, from this recoupling to the computation 



of the coefficient functions can be reduced to a Mathematica program. 
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